AFIT/GE/EE/78 -39 


<0 

(M 

O 


Q_ 

O 

C_> 


SIMULTANEOUS  NON-HOMOGENEOUS 
MULTIPLE  ACCESS  TECHNIQUES 


THESIS 


AFIT/GE/EE/78-39  Pedro  Rodriguez 

Capt.  USAF 


Of/ 


Approved  for  public  release;  distribution  unlimited. 

V"’ 

- or.  H ■ HT.  ?S , KaJ  or , Um5F 

Director  of  Info. nation  V_)  / i Pf 


V 


CO 


7T- 


AFIT/GE/EE/78 - 39 


SIMULTANEOUS  NON -HOMOGENEOUS  MULTIPLE 


ACCESS  TECHNIQUES. 


™ESIS. 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  Training  Command 
in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


Ci Pedr<^kodr igue z J B.S. 


Capt . 


USAF 


Graduate  Electrical  Engineering 


Dec 


Approved  for  public  release;  distribution  unlimited. 


cl&- 


Preface 


I 


The  purpose  of  this  thesis  was  twofold.  The  first  was  to 
satisfy  the  School  of  Engineering  thesis  requirement.  The 
second  was  to  investigate  a satellite  transponder  system  which 
is  already  in  the  development  stage  to  determine  the  best  use 
of  that  system  when  accessed  by  the  equipment  presently  in  use 
and  planned  to  become  operational  in  the  near  future.  A study 
of  a communications  control  system  was  originally  part  of  the 
topic.  Control  systems  are  presently  being  investigated 
feverishly.  Most  articles  which  have  been  published  recently 
conclude  that  the  control  function  should  be  placed  at  the 
satellite. 

I am  extremely  grateful  to  Capt.  Thourot  from  Hq.  AFCS, 
who  submitted  a proposal  for  consideration  as  a thesis  topic 
from  which  this  thesis  topic  was  extracted.  I wish  to  thank 
my  readers,  Maj . Cappinella  and  Capt.  Hadley.  I wish  to  ex- 
press appreciation  to  Asst.  Prof.  Joseph  Carl  of  the  Department 
of  Electrical  Engineering,  for  his  helpful  guidance,  and  coun- 
seling during  this  study.  I owe  a special  thank  you  to  Mrs. 
Weber,  who  typed  the  final  copy. 

I wish  to  express  fatherly  love  to  my  daughters  who  often 
experienced  my  wrath  brought  about  by  school  stresses,  and 
brought  about,  they  must  have  felt,  for  no  apparent  reason. 

I wish  to  thank  my  wife  for  her  patience  and  consideration 
without  which  this  work  would  never  have  been  completed.  I 
dedicate  this  effort  to  my  soon- to-be-born  child. 
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Abstract 

The  possibility  of  using  simultaneous  non-homogeneous 
multiple  access  techniques  through  a common  transponder 
channel  was  investigated.  Three  configurations  were  analyzed: 
(1)  a single  time  division  multiple  access,  (2)  several  fre- 
quency division  multiple  accesses,  and  (3)  the  combined 
accesses  of  configurations  (1)  and  (2) . Each  of  the  configura- 
tions was  studied  in  three  disjoint  segments:  Cl)  the 
bandwidth  and  shape  of  the  power  spectrum  for  various 
modulating  waveforms,  (2)  the  probability  for  error  for 
linear  PSK  and  non-linear  PSK  systems,  and  (3)  the  distortion 
effects  that  a signal  experiences  through  a non-linearity. 

Each  segment  was  studied  in  detail  independently  of  the  others. 
The  comprehensive  analysis  accomplished  on  each  segment,  when 
joined  into  a system,  allows  the  end-to-end  channel  characteri- 
zation.^A  link  calculation  was  accomplished  for  a single  TDMA 
access,  because  it  was  found  to  be  the  optimum  system  based  on 
the  analysis  of  the  segments.  Also,  the  results  obtained 
indicate  that  the  benefits  derived  by  using  non-homogeneous 
accesses  through  a common  transponder  channel  are  far  less 
than  the  total  degradation  the  system  experiences.  A recom- 
mended alternative  is  provided. 
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SIMULTANEOUS  NON -HOMOGENEOUS  MULTIPLE 


ACCESS  TECHNIQUES 

I . Introduction 

Background 

The  Defense  Satellite  Communications  System  (DSCS)  Program 
currently  includes  research  and  development  of  and  plans  to 
phase-in  time  division  multiple  access  (TDMAj  equipment.  It 
is  anticipated  by  Headquarters  Air  Force  Communications 
Service  (AFCS)  that  during  the  middle  1980's  a transition  will 
occur  in  the  use  of  satellite  multiple  access  techniques  from 
frequency  devision  multiple  access  (FDMA)  to  TDMA.  The  use 
of  FDMA  will  not  immediately  cease  when  this  transition  takes 
place,  but  will  continue  for  the  duration  of  the  FDMA  equip- 
ment's useful  life.  It  is  for  that  reason  that  this  study  was 
accomplished  to  determine  possible  problems  which  might  be 
encountered  when  both  multiple  access  techniques  are  used  simul- 
taneously through  the  same  satellite  transponder. 

Problem 

The  problem  that  was  addressed  was  one  of  intermodulation 
products  and  signal  suppression  factors.  For  TDMA,  a trans- 
ponder is  operated  at  saturation  to  take  full  advantage  of  the 
traveling-wave-tube  (TWT)  output  power.  For  FDMA,  since 
several  carrier  frequencies  are  present  at  the  transponder 
input,  the  power  level  of  each  carrier  must  be  monitored  and 
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controlled  to  minimize  intermodulation;  nevertheless,  inter- 
modulation products  still  result  and  reduce  the  available 
output  power  of  the  desired  signals.  When  FDMA  and  TDMA  are 
used  simultaneously  through  the  same  transponder,  the  weak- 
nesses inherent  in  each  are  present.  The  purposes  of  this 
study  are  to  investigate  those  weaknesses  to  find  an  alterna- 
tive or  alternatives  that  make  the  best  use  of  the  transponder's 
available  power  and  bandwidth,  and  to  provide  recommendations 
as  to  the  best  alternative. 

Scope 

This  study  dealt  almost  entirely  with  the  analysis  of  a 
transponder  accessed  simultaneously  by  FDMA  and  TDMA  techniques 
to  determine  what  degradations,  if  any,  occur  in  the  system. 
Theoretical  concepts  that  are  closely  related  to  the  problem 
under  study  are  included  in  this  report.  This  study  does  not 
address  the  communications  control  or  synchronization  timing 
problems,  since  those  problems  will  be  solved  during  the  design 
of  the  TDMA  equipment  subsequent  to  this  study.  The  DSCS  III 
space  segment  and  the  SHF  Ground  Mobile  Forces  (GMF)  terminals 
are  included  as  constraints  on  the  analysis,  but  in  reality, 
the  analysis  is  equally  applicable  to  other  similar  systems. 

Assumptions 

Three  assumptions  were  made  and  are  stated  as  follows: 
the  noise  which  corrupts  the  transmission  channel  is  white 
Gaussian  noise  which  has  zero  mean  and  has  a flat  power  spec- 
trum, the  noise  bandwidth  is  much  smaller  than  the  center 
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frequency  of  the  bandpass  filter,  and  the  bandpass  filter 
which  precedes  and  the  zonal  bandpass  filter  which  follows 
the  limiter  are  wide  enough  to  pass  the  desired  signal  with 
negligible  distortion. 

Summary  of  Current  Knowledge 

Of  the  three  generic  satellite  multiple  access  techniques, 
TDMA,  FDMA,  and  spread- spectrum  multiple  access  (SSMA) , TDMA 
achieves  the  best  power  and  bandwidth  efficiency.  The  power 
efficiency  advantage  of  TDMA  results  from  the  lack  of  inter- 
modulation distortion,  since  only  one  carrier  is  present  at 
the  transponder.  Overhead  bits  which  are  required  for  guard 
time,  synchronization,  and  network  control  functions  reduce 
the  power  efficiency.  In  FDMA,  power  losses  result  from 
intermodulation,  which  is  created  by  multiple  accesses  occupy- 
ing the  transponder  simultaneously,  and  suppression  of  one 
signal  by  another.  To  reduce  intermodulation  the  transponder 
is  backed-off  from  saturation  to  a nearly  linear  operating 
region,  but  this  results  in  power  efficiency  loss  and  capacity 
loss  due  to  reduced  down  link  effective  radiated  power  (ERP) . 
Another  FDMA  loss  occurs  in  the  use  of  real  filters  which  have 
non-ideal  roll-off  characteristics,  so  that  a guard  band  has  to 
be  provided  between  accesses  to  keep  the  losses  due  to  cross- 
talk to  a tolerable  level. 

Approach 

The  problem  was  analyzed  in  three  configurations  which  are 
as  follows:  1)  a single  TDMA  carrier  and  its  optimum 
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performance,  2)  several  FDMA  carriers  and  their  optimum 
performance,  and  3)  a compromise  configuration  consisting  of 
a single  TDMA  signal  and  several  FDMA  signals.  The  perform- 
ance parameters  were  found  for  each  of  the  three  cases  listed 
above.  However,  before  the  performance  parameters  could  be 
determined  it  was  necessary  to  segment  the  problem  into 
disjoint  parts  and  to  analyze  those  parts  in  detail. 

The  individual  parts  of  the  problem  as  finally  segmented 
follow:  1)  the  derivation  of  the  power  spectrum  for  phase- 

shift-keyed  modulation,  2)  the  analysis  of  a bandpass  limiter 
when  it  is  in  a hard  and  soft  configuration,  and  3)  the 
detection  of  the  PSK  signal  through  linear  and  ron-linear 
limiters.  Once  the  individual  parts  were  analyzed,  the  pieces 
were  joined  and  studied  as  a whole.  Link  calculations  were 
then  possible  since  the  parameters  for  each  segment  were  known. 

In  addition,  the  several  cases  of  concern  were  analyzed  to 
determine  the  feasibility  of  using  simultaneous  non-homogeneous 
multiple  access  techniques  through  a common  transponder  channel. 

A large  amount  of  information  is  available  on  this  subject. 

Many  of  the  figures  included  herein  were  generated  from  equa- 
tions found  in  the  literature  and  were  numerically  evaluated 
by  the  computer.  The  results  obtained  were  then  plotted  by 
the  calcomp  plotter.  Most  of  the  equations  found  in  the 
literature  involved  Bessel  functions,  so  graphical  results  are 
provided  as  visual  aids  to  make  clear  th>:  significance  of  those 
equations . 
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II.  Power  Spectrum 


The  power  spectrum  of  digital  angle  modulated  signals  can 
be  calculated  in  many  ways.  The  direct  way  of  calculating 
the  power  spectrum  is  to  first  find  the  Fourier  transform  of 
a sample  of  the  signal  on  a finite  interval.  This  approach 
was  used  here.  However,  the  double  Fourier  transform  had  to 
be  used  since  the  signal  under  consideration  was  found  to  be 
non-stationary . 

Digital  binary  phase-modulated  signals  can  be  represented 
as 

N 

x(t)  = I b cos(w  t)p(t-nT)  (1) 

n=0  n c 

where  the  bn's  are  identically  independently  distributed  ran- 
dom variables  which  each  take  values  +1  and  -1  with  equal 
probability,  and  p(t)  is  the  binary  pulse.  In  this  analysis 
p(t)  is  a rectangular  pulse  of  time  duration  T. 

Derivation 

The  power  spectral  density  of  the  random  process  is  the 
Fourier  transform  of  the  autocorrelation  function  which  is 
obtained  from 

Rx(tl»  t2)  Sx(f)  (2) 

where 

W V ■ (5) 
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Also,  autocorrelation  is  equal  to  the  expected  value  of  the 
random  process  sampled  at  two  arbitrary  time  instances.  It 
is  expressed  as 


Rx(tl»t2^  = E[x(t1)x(t2)  ] 


/ N 

Rx^tl»t2^  = El  E0bn  cos  Vl  P(ti"nT) 

N 

• E b cos  w t~  p(t?-mT) 
m=0  m C-  l l 

N N 

= cos  w t.  cos  w t-E  E E b b • p (t. -nT) p (t7 -mT)  (4) 

c 1 c A=0  m=0  n m 1 £ 


But 

E t t>_b_  ] = \?'  ™ ^ =6  , the  Kronecker  delta,  so 

1 n m J ( 1 , m = n mn’ 

N 

Rx^tl,t2^  = cos  wctl  cos  wct2  E P(ti"nT^)  * P C 1 2 "nT)  (5) 


Now,  since  p(t^-nT)  and  p(t2~nT)  are  both  sample  values 

of  the  same  shifted  rectangular  pulse  (which  is  non-zero  only 

on  an  interval  of  length  T) , their  product  is  non-zero  only 

if  t^  and  t2  lie  in  the  same  signaling  interval  of  length  T 

for  exactly  one  value  of  n.  In  that  case,  and  that  case  only, 

the  product  has  the  same  value  as  the  square  of  the  pulse 
2 

amplitude,  say  A ; otherwise  the  product  is  zero.  Thus, 


P-y  (t-ljtO  = 


A cos  w t,  cos  w t,,  if  t.  and  t0  lie  in 
cl  c 2 1 Z 

the  same  signalling  interval 
0,  else 
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w 

A process  is  stationary  in  the  wide  sense  if  its  expected 
value  is  a constant  and  its  autocorrelation  depends  only  on 
the  difference  x = t^-t2.  TEe  Process  in  Ecl  (5)  has  values 
other  than  zero  when  the  selected  arbitrary  time  instances 
fall  in  the  same  signaling  interval  and  is  identically  zero 
when  the  selected  arbitrary  time  instances  fall  in  different 
signaling  intervals.  Since  the  cosine  terms  are  also 
expressible  as 

cos  w(.t1  cos  wct2  = y[cos  wc(ti  + t2^  + cos  wc(-tl't2^ 

It  is  obvious  that  the  random  process  under  consideration  does 
not  have  an  autocorrelation  function  which  is  solely  a function 
of  t^-t2.  Hence,  the  process  is  not  stationary,  even  in  the 
wide  sense.  Its  spectrum  cannot  be  obtained  from  the  Fourier 
transform  of  the  autocorrelation  function. 

Because  the  autocorrelation  functions  of  non-stationary 
processes  are  functions  of  two  variables,  in  order  to  apply 
transform  techniques,  two-dimensional  Fourier  transforms  are 
needed.  The  double  Fourier  transform,  r(w^,w2),  of  Rx(t^,t2) 
can  also  be  defined.  The  double  Fourier  transform  can  be 
represented  as  (Ref  12:466) 

r(w1,w2)  = E[X(w1)X*(w2)1  (7a) 

- f?t?  Rx(t1,t2)e'j(wltrw2t2)dt1dt2  (7b) 

Thus,  since  the  process  at  hand  is  non-stationary,  this  more 
general  approach  for  finding  the  power  spectrum  must  be  used. 
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The  random  process  will  be  first  transformed  into  the  fre- 
quency domain  using  the  Fourier  transform.  The  expected  value 
of  the  new  random  process  multiplied  by  its  conjugate  yields 
the  expression  that  represents  the  power  spectrum  of  the 
random  process.  The  random  process  in  the  frequency  domain 
is  given  by 

X(f)  = ? bn  i{6(f-fc)  + 6(f+f  )}  * [P(f)-ej27rfnT]  (8a) 

n=0 

N 

X(f)  = Z b,  i[P(f-f  ) + P(f+f  )]ej2TrfnT  (8b) 


and  its  conjugate 


X (f ) = Z bm  £[P*(f-fc)  + P* (f+fc) ] e 
m=  0 


■ j 2irfmT 


Multiplying  X(f)  by  X*(f),  taking  the  expected  value  of  the 
product,  and  assuming  positive  and  negative  frequency  terms 


do  net  overlap  yields 


N N 


E[X(f)X*(f)l  = T 2 2 [ | P(f-f _) | 2 + |P(f+f.)|2] 

4 n=0  m=0  c 


• E [b  b ] .ej2Tt£T(n_m^ 
n mJ 


But  E[b  b 1 = 6 , so  that  if  the  cross-product  terms  like 

L n mJ  mn 

P(f-f  )P*(f+f  )--it  can  be  argued  that  f is  large  enough 
c c c 

that  positive  and  negative  frequency  terms  do  not  overlap--, 


then 


E[X(f)X* (f) ] = h § |P(f-fc)|2  ♦ |P(f+fc)|2]  (10a) 
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Fig  1.  Typical  Rectangular  Modulating  Waveform 
or 

E[|x(f)|2]  = 5r[|P(f-fc)  |2  + |p(f+fc)|2]  (iob) 

The  assumed  form  of  the  modulating  signal  is  shown  in  Fig  1. 

The  transform  of  a rectangular  pulse  is  easily  obtained  from 

a transform-pair  table.  Most  often  the  rectangular  pulse  is 

found  shifted  T/2  to  the  left  in  the  transform  table.  However, 

since  the  only  term  sought  here  is  the  magnitude,  the  e^  and 
- i 4* 

e J terms  will  cancel  each  other.  The  magnitude  of  the 
Fourier  transform  of  p(t)  can  be  found  from 

p(t-y)  P(f)  = sin  irfT  £ AT  sinc( fT)  (11) 

| P(f ) | 2 = |P(f)|2  = A2T2sinc2(fT)  (12) 

hence , 

E [X(f ) X* (f) ] = NA2T2/ 4 [sinc2T(f-fc)+sinc2T(f+fc)]  (13) 
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Comparison  of  Spectrum 

The  power  spectrum,  as  derived  and  given  by  Eq  (13) , can 
be  related  to  the  results  that  Glance  (Ref  5:2869)  derived  for 
a rectangular  pulse  of  duration  t<T.  His  results  are  given 
below.  The  continuous  part  of  the  spectrum  is  of  the  form 

G(f ) = (T/4)(t/T)2  {sin[2TT(fc-f)  (x/2)]/(fc-f)  (ttt)}2  (14) 

but  here  t=T.  With  minor  mathematical  manipulations,  Eq  (14) 
can  be  written  as 

G (f ) = (T/4)  (sin[2TT(fc-f)(T/2)]/(fc-f)(TTT)}2  (15) 

which  is  identical  to  Eq  (13) . 

Having  pointed  out  that  the  results  derived  here  for 
BPSK  modulation  agree  with  those  of  Glance,  the  results  found 
in  Glance's  article  are  used  hereafter  since  they  are  in  a 
more  general  form.  The  trapezoidal  waveform  was  used  since  it 
better  represents  a real  baseband  modulating  signal.  Now  the 
modulating  waveform  appears  as  in  Fig  2.  It  is  composed  of 
several  signals  that  have  been  time  division  multiplexed  into 
a continuous  baseband  data  stream. 

Binary  Phase-Shift  Keyed 

Glance  (Ref  5:2870-2871)  has  calculated  the  power  spectrum 
for  antipodal  PSK  modulation  by  a trapezoidal  waveform.  The 
spectrum  has  a continuous  part  and  a discrete  part  as  follows: 
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n2 


Gc<f)  " T<T> 


t/T  sin  [ (fc-f)  (ttt)  1 , S rrc  n. 
*-*-  + y cos[(£c-£)  (ttT)] 

[ C£c-£)  (S/2)  ] 2 - (tt/4)  2 


(16) 


for  the  continuous  part,  and 

Gd(£)  " T(f2)(r)2  2 — q.pg2i2Trm  S/T)  & (f  +f •*•!«•)  a?) 

41  m=-oo[(irm^)z  - (tt/4)Z]Z  C 1 


for  the  discrete  part.  Because  this  representation  makes  x(t) 
have  zero  mean,  there  is  no  discrete  part  to  its  spectrum. 

The  power  spectrum  given  by  Eq  (16)  decreases  as  1/f4  as  £-*•«». 
The  spectral  density  and  modulating  signal  representation  for 
the  continuous  part  are  shown  in  Fig  3.  The  power  spectra 
has  a lower  amplitude  and  greater  bandwidth  as  the  modulating 
rate  is  increased.  This  will  become  a very  important  relation- 
ship later  when  the  power  in  the  carrier  is  related  to  the 
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Fig  3.  Spectral  Density  vs  Frequency,  (a)  Phase 
Modulation  - 10  Mb/s  (b)  Phase  Modulation- 

5 Mb/s 
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energy-per-bit-E^ . The  energy-per-bit  is  defined  as 


E,  = A2T/2  (18) 

D 

where  A is  the  amplitude  of  the  sinusoidal  constant  envelope 
r-f  carrier  and  T is  the  information  bit  time  duration.  For 
a BPSK  signal  with  phase  modulation  rate  of  10  Mb/s,  the 
required  bandwidth  is  20  MHz.  The  bandwidth  requirement  is 
linearly  related  to  the  modulation  rate.  For  BPSK,  the 
bandwidth  requirement  is  twice  the  modulation  rate.  For  QPSK, 
the  bandwidth  is  related  to  the  modulation  rate  on  a one  to 
one  basis. 
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III . Bandpass  Limiter 

There  are  two  effects  that  a signal  passed  through  a 
limiter,  such  as  a TWT,  experiences.  These  effects  are  signal 
suppression  and  distortion.  For  instance,  the  result?  of  this 
chapter  show  that  the  suppression  a signal  experiences  in  the 
presence  of  n signals  which  have  equal  magnitude  is  -1.05  dB. 
Additionally,  their  interaction  in  the  non-linear  device  creates 
intermodulation  products  (IMP's)  which  further  degrade  the 
desired  signal  terms.  How  much  signal  suppression  and  dis- 
tortion a non-linearity  introduces  into  a communications 
channel  is  investigated  in  this  chapter. 

Moreover,  a method  is  presented  to  reduce  the  harmful 
effects  introduced  by  the  non-linearity.  Also,  a method  is 
presented  to  control  the  harmful  effects  of  IMP's.  These 
methods  have  been  used  extensively  in  the  past.  They  use  both 
the  careful  control  of  the  drive  levels  to  the  input  of  the 
TWT  and  judicious  frequency  assignment  of  the  r-f  carriers. 

Hard-Limiter  Suppression 

An  analysis  of  the  output  signal  from  the  bandpass  limiter 
shows  the  relationship  between  the  input  and  output  signal-to- 
noise  ratio.  The  input-output  relationship  is  determined  for 
coherent  and  incoherent  demodulation.  The  results  obtained 
for  each  of  the  two  cases  are  compared. 


14 


Communications  Path  Model.  Figure  4 is  a block  diagram 
showing  the  bandpass  limiter.  In  that  figure,  h^  is  a band- 
pass filter  which  has  center  frequency  wc , while  h2  is  a 
zonal  filter  which  selects  the  first  zone  components  only. 

It  is  assumed  that  the  bandwidth  of  the  filter  h^  is  much 
less  than  its  center  frequency. 

The  PSK  signal  plus  noise  is  defined  to  be  of  the  form 

x(t)  = s ( t ) + n(t)  (19) 

where 

s ( t ) = A cos  [w£t  + <J>(t)]  (2C) 

where  4>(t)  is  a random  phase  process  which  takes  values  of  0 
and  u,  and  is  the  information  bearing  signal.  The  relation- 
ship to  the  notation  in  Chapter  II  should  be  clear.  The 
noise  is  assumed  to  be  white  Gaussian  noise  expressed  as 

n(t)  = nc(t)  cos[wct  + <{>(t)]  - ns(t)  sin  [wct+<f>  (t)  ] (21) 

where 


B(ns(t)]  - E[nc(t)]  - 0 

(22) 

E tns (t) ] - E[n|(t)]  - ol/2 

(23) 

°n  " NoBi 

(24) 

is  the  input  bandwidth  of  filter  h^  in  Fig  4 . The  ideal 
limiter  function  is  defined  as 
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Fig  4.  PSK  System  Block  Diagram  Model 

l+lt  x(t)  > 0 

yT  (t)  = sgn  [x  (t ) ] = < (25) 

l-l,  x (t)  < 0 

Since 

E[x(t)]  = A E{cos  [v\'ct+(j)  (t)  ] } (26) 

the  envelope  of  y^(t)  = A.  Using  trigonometric  identities, 
the  input  process  x(t)  can  be  written  equivalently  as 

x (t)  = a (t)  cos  [wct+<J>  (t)+y(t)  ] (27) 

where 

a (t)  = ( [A+nc (t) ] 2 + ng (t) } 1/2  (28) 

and 

y(t)  = tan'1{ns(t)/[A+nc(t)] } (29) 

The  signal  given  by  Eq  (27)  has  been  analyzed  by  several 
individuals.  Middleton  (Ref  11:144)  established  the  input- 
output  time  function  relationship  for  the  limiter  using  the 
complex  Fourier  transform.  The  output  from  the  limiter  as 
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derived  by  Middleton  is  given  by 


yL(t)  = l/2u  fQ  f(j£)  exp  [j£x(t)]  d £ (30a) 

where  f(j£)  = 2/j  is  the  Fourier  transform  of  the  dynamic 
path  represented  by  the  limiter.  By  expanding  the  exponen- 
tial function  in  a series,  the  limiter  output  can  be  written 
as  (Ref  16:43) 

yL(t)  = 4/tt  {cos[wct+Y(t)+<Kt)] 

+ 1/3  cos  [3wct  + 3y  (t)  + 3<{>  (t)  ] + ...}  (30b) 

The  output  of  the  limiter  is  represented  by  Eq  (30b)  as 
a series  of  phase-modulated  sinusoids.  This  series  begins 
with  frequency  wc  and  extends  to  all  odd  multiples  of  wc> 
Since  h2  is  a zonal  filter  which  selects  the  first  zone  com- 
ponents only,  the  first  zone  output  is 

yx(t)  A 4/tt  cos  [wct+y  (t)  +<{>  (t)  ] (31) 

In  reality,  a physical  bandpass  filter  can  only  approximate 
the  required  output.  However,  because  of  the  narrowband 
restriction  which  was  stated  earlier  in  this  section,  the 
error  is  small.  The  expected  value  of  the  power  in  the  first 
zone  components  is  (Ref  16:43) 

E(yx  (t)  2}  - 8/tt2  (32) 

Incoherent  Detection.  The  system  under  consideration  is 
assumed  to  be  coherent.  However,  the  incoherent  results  are 
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included  here  for  completeness.  These  results  might  prove 
useful  later  on  or  will  serve  to  illustrate  a point  when 
compared  to  the  coherent  results. 


Using  Davenport's  method  (Ref  3)  which  amount  to  evaluating 
the  amount  of  carrier  component  through  formation  of  the 
second  moment  function  and  averaging  over  the  distributions 
of  <(>(t)  and  n(t),  Springett  and  Simon  (Ref  16)  derive  the 
bandpass  limiter  (BPL)  input-output  relationship.  Their 
results  follow  in  the  remainder  of  this  section.  The  first 
zone  output  is 

y1(t)  = u(p)  cos  [wct+(f>  (t)  +n ' (t)  ] (33) 

where  n' (t)  consists  of  the  interaction  of  signal  and  noise 
(sxn) , and  noise  and  noise  (nxn) , and  u(p)  is  defined  below. 

The  limiter  input  SNR  is  p,  defined  as 

(S/N)i  A p = A2/2NQBi  (34) 

The  carrier  component  is  now  represented  as  a function  of 

P by 

2 

u ( p)  = (2p/7T) 1/2  e'p  {I0(p/2)  + I1  (p/2)  } (35) 

where  I and  I,  are  the  modified  Bessel  functions  of  the  first 
kind  of  order  zero  and  one,  respectively.  The  noise  ana 
signal  power  are  now  given  as  functions  of  u(p),  which  is  a 
function  of  a(p).  The  function  a(p)  is  commonly  referred  to 
as  the  signal  component  suppression  factor.  The  output  noise 
power  is  given  by 
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(36) 


Pn  = 8/u2-u(p) 2 = 8/n2 [l-a(p) 2] 
and  signal  power  is  given  by 

Ps  = u(p)  2 = 8 / tt 2 a (p)  2 (37) 

which  give  an  output  SNR  as 

2 

(S/N)*  = a (-P-3-  -j  ■ (38) 

° [l-a(p)  ] 

The  superscript  I denotes  ’'incoherent,"  and  the  subscript  0 
denotes  "output." 

Coherent  Detection.  The  reference  signal  used  for  coher- 
ent detection  is  assumed  to  be  of  the  form 

r (t , 6)  = 2 cos[wct  + e(t)+4>(t)]  (39) 

where  9 takes  values  between  (tt/2  , -tt/2)  and  is  the  demodula- 
tion reference  angle.  The  phase  angle  cf>  ( t ) is  taken  as  known 
a priori.  The  signal  and  noise  are  as  previously  defined  in 
Eqs  (20)  and  (21)  from  which  the  output  y^(t)  of  a BPL  is 
derived  as 

yx(t)  = 4/ it  cos[wct+Y'  (t)]  (40) 

where 

A sin  <p  (t)  +n  (t) 

Y’(t)  = tan  } (41) 

A cos  <t>  (t)  +nc  (t) 

After  substituting  the  proper  expressions  for  sin  Y'(t)  and 
cos  Y’(t),  the  demodulated  output  obtained  is  written  as 


19 


(42) 


Z(9)  = 2- 


A + nc  ns 

2 — 2 1/2  cos  0 + 2 — 2” ITT  sin  e 

[ (A+nc) 2+ng] 17  ^ [(A+nc)Z+n;]1/Z 


The  first  two  moments  of  Eq  (42)  are  now  required  to  evaluate 
the  SNR.  The  first  moment  uz  is  defined  as 


uz  = E{ z (0 ) } = z fz  dz 


(43) 


while  the  second  moment  is  defined  as 


E{z(0)2}  = Cjz1  fz  dz 


(44) 


These  two  moments  are  then  used  to  express  the  variance  as 

(45) 


y ? y 

oj  = E[zZ]  - E [z] ~ 


Because  it  has  been  assumed  that  the  noise  is  white  Gaussian 
noise,  the  evaluation  of  the  expectations  yield  (Ref  16:43-45) 

,1/2  -P/2 


uz  = ( 2p/tt) 


e [I0(p/2)+I1  CP/2)  ] cos  0 (46) 


and 


r2  = 8 / tt 2 [cos2  0 - (l-e”p/2p)  cos  20] 


u. 


(47) 


Since  the  input  and  output  SNR  are  the  quantities  of  interest, 
these  SNR's  are  given  by 


S 0 uz 

= T 

°z 


(48) 


at  the  output  of  the  BPL . The  output- input  ratio  is 
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Fig  5.  -Bandpass  Limiter  Ratios  of  SNRs 


$1 

2(ff}i 


(49) 


The  bandwidth  change  of  1/2  is  required  due  to  the  bandpass 
to  lowpass  transformation.  The  input  SNR  is  p.  Finally,  the 
SNR  for  the  coherent  case  is  a function  of  p when  6 = 0 is 
expressed  as 

CS/Nlg  _ j a -p[I0tP/2)  + I^P/2)]2 (50) 

2(S/N)i  2 (4/Ti)(l-[(l-e'p)/2  p]}-pe"p[Io(p/2+I1(p/2)]Z 

% 

and  when  6 = u/2  is  expressed  as 
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(51) 


(S/N)^/2  _ 1 e 'P[I0(p/2)+I1(p/2)]2 
2(S/N).  1 2/tt  [(l-e  ~p)/p  ] 

In  Fig  5,  the  relationship  between  the  input-output  ratio 
for  9 = ff/2  is  graphically  depected.  Included  in  this 
section  are  the  incoherent  and  coherent  SNR  relationships. 

The  complete  end-to-end  hard-limited  channel  for  a TDMA  or 
an  FDMA  access  with  PSK  modulation  can  be  characterized  by 
parameters  such  as  P , SNR,  BW,  and  bit  rate. 

Soft-Limiter  Suppression 

The  analysis  of  the  soft-limiter  follows  that  of  the 
hard-limiter.  Springett  and  Simon  (Ref  16)  have  found  that 
when  a hard-limiter  is  followed  by  a coherent  detector,  the 
strong  signal  asymptotic  values  no  longer  agree  with  the 
incoherent  results  derived  by  Davenport.  Those  results  were 
found  to  apply  equally  well  to  a coherent  soft-limiter. 

The  signal  and  noise,  as  previously  defined  in  the  pre- 
ceding section  and  given  by  Eq  (19),  are  once  again  the  input. 
The  model  for  the  soft-limiter  is  given  in  Fig  6 and  replaces 
the  ideal  limiter  in  Fig  4 given  by  Eq  (25).  The  soft-limiter 
representation  seen  in  the  figure  is  that  for  the  error 
function  representation. 

The  narrowband  waveform  which  is  represented  by  Eq  (20), 
when  fed  to  a memoryless  non-linear  device,  yields  an  output 
v(t),  which  is  a function  of  the  input.  Before  proceeding 
further,  the  input  was  written  in  a more  convenient  form  as 
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Fig  6.  Soft-Limiter  Model  (Represented  by 
the  Error  Function 


x(t)  = a (t)  cos[wct  + <f>  + y (t)  ] (52) 


The  input  signal  and  noise  are  identical  to  those  given  by 
Eqs  (27),  (28),  and  (29).  Letting  a = wct  + 4>  + Y(t),  the 
input  is  now  expressed  as 


x(t)  = a(t)  cos  a 


(53) 


Since  the  output  is  an  even  periodic  function  cf  a,  it  is 
represented  by  the  Fourier  series  (Ref  2:398-399) 


v(a  cos  a)  = 1/2  vQ(a)  + v^a)  cos  a + . . . (54) 

The  expansion  is  valid  regardless  of  how  a(t)  and  a vary  in 
time.  The  vm(a)  terms  of  Eq  (54)  are  given  by  the  Chebyshev 
transform.  These  terms  are  called  the  m-th  order  of  the 
transform.  The  first  order  Chebyshev  transform,  v, (a)  is  the 
only  frequency  component  that  is  passed  by  the  zonal  bandpass 


23 


filter  which  follows  the  non-linearity.  The  variable  a(t) 
is  non-negative,  and  is  defined  by  the  expression  of  Eq  (28). 
The  non-linearity  for  a soft-limiter  is  expressed  as 


y(t) 


Ktt 


1/2 


L er£  CiirxrtTi 


(55) 


where  L is  the  peak  limiter  output  and  K is  the  slope  of  the 
transfer  characteristic  at  x = 0.  Hence,  the  output  of  the 
soft-limiter  as  determined  by  Blachman  (Ref  8:399)  is  given 
by  the  expression 

-a2 

1 a 2 2 

7l(t)  = L(f)7  ae"1”  {I^)  + I^)}  (56) 

The  independent  variable  t has  been  suppressed  above.  The 
expression  of  Eq  (56)  was  numerically  evaluated,  and  those 
results  are  provided  in  graphical  form  in  Fig  7.  The  ampli- 
tude of  the  signal  approaches  its  maximum  value  (limit)  for 
SNR>4.  At  first  glance,  it  appears  that  the  signal  exceeds 
unity,  and  hence,  the  output  appears  to  be  amplified. 

However,  that  is  not  the  case.  The  amplitude  obtained  at 
the  output  of  a hard-limiter  must  be  multiplied  by  the  sup- 
pression factor  tt/4.  Thus,  the  amplitude  reduces  to  a 
maximum  output  value  of  unity.  Consequently,  v^(a)  describes 
the  envelope  in  the  fundamental  zone.  The  envelope  was 
determined  from  the  results  of  Blachman  (Ref  2)  who  made 
use  of  the  Chebyshev  transform  to  derive  his  results. 

A coefficient  must  now  be  introduced,  since  it  was  ignored 


completely  in  the  previous  steps.  Because  the  error  function 
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Fig  7.  Envelope  Magnitude  at  the  Output  of  the 
Soft-Limiter  versus  Input 
Signal -to-Noise  Ratio 

1/2 

was  defined  with  the  argument,  [K  /2L  x(t)],  the  coeffici- 
ent that  precedes  x(t)  must  be  accounted  for.  It  is  obvious 
that  a(t) , as  previously  defined,  must  be  modified  to  the 
expression 

\ i 

aft)  - { [A  ♦ nc(t)l2  ♦ n2  ft) } 2 (57) 

which  simply  says  that  the  peak  value  of  the  limiter  and 
slope  at  x = 0 are  variables.  They  are  defined  specifically 
for  a specific  T'.'/T . 

The  output  y^(t),  which  is  similar  to  that  of  Eq  (40),  is 
then  multiplied  by  the  coherent  reference,  r(t)  = cos(wct  + 6), 
and  low-pass  filtered.  After  the  expressions  for  sin  y(t)  and 
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cos  y (t)  are  substituted  into  the  resultant  output,  the  output 
is  identical  to  that  of  Eq  (42),  but  it  must  be  multiplied  by 
the  coefficient  -K1/,2/2L. 

The  first  two  moments  of  the  output  are  required  to 
evaluate  the  SNR.  The  first  moment  uz  is  defined  in  Eq  (43), 
while  the  second  moment  is  defined  in  Eq  (44)  . From  these 
two  moments,  the  variance  was  found  using  Eq  (45).  The  process 
used  to  derive  the  first  and  second  moments  is  not  summarized 
here.  Instead  the  reader  is  asked  to  read  (Ref  9)  if  he 
wants  that  information.  Proceeding  further  and  using  Lesh's 
(Ref  9)  results,  the  input-output  signal-to-noise  ratios  for 
a coherent,  phase  detector  is  given  by 


snrJ/2 

2SNRi 


"uz | 9 =0 

R ' { E [ z 2 ] - u2} 


(58) 


where  0 = ir/2.  R'  is  the  SNR  at  the  input  after  projection 
onto  the  coherent  reference  axis  and  is  given  by 

R'  = 2R  cos2  9 (59) 


which  further  reduces  the  expression  of  Eq  (58)  to  an  expres- 
sion where  R1  is  replaced  by  2R. 

The  behavior  of  a phase  detector  with  detection  angle, 

0 = tt/2,  is  expressed  in  terms  of  input-output  SNR  relation- 
ships in  Fig  8.  The  parameter  B = L/KA  in  that  figure  can  be 
interpreted  as  the  ratio  of  the  peak  limiter  output  to  the 
peak  signal  output,  if  the  limiter  is  a linear  amplifier  with 
gain  K.  When  B = 0,  the  error  function  represents  a 
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Fig  8.  Soft  Bandpass  Limiter  Ratios  of  SNRs 


hard-limiter.  The  curves  in  the  figure  were  extracted  from 
another  figure  found  in  ( Ref  9).  The  curve  for  8 = 0 is 
the  actual  curve  that  is  generated  for  a hard-limiter  as  a 
function  of  SNRs.  The  phenomenon  that  is  depicted  by  the 
three  curves  is  obvious.  As  0 gets  larger,  the  system  behaves 
as  a linear  amplifier.  As  3 gets  smaller,  the  system  behaves 
as  a soft-limiter.  The  system  becomes  a hard-limiter  when 
0 = 0. 


Spectrum  Distortion 

The  problem  under  consideration  is  the  amplification  of 
several  signals  in  a common  device  in  which  limiting  occurs. 
This  problem  is  of  particular  importance  in  satellite  communi- 
cation systems  when  more  than  one  access  is  required.  The 
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method  which  was  derived  from  Rice  (Ref  3)  and  used  here  in 
determining  the  output  of  a limiter  follows  that  of  Davenport 
and  Root  (Ref  4:277-311).  It  consists  of  computing  the  auto- 
correlation function  at  the  limiter  output  and  then  taking 
a Fourier  transform  to  obtain  the  spectrum.  The  derivation 
of  the  autocorrelation  and  spectrum  is  summarized  and  follows. 
The  Wiener-Khintchine  theorem  states  that  the  Fourier  trans- 
form of  the  spectral  density  yields  the  autocorrelation 
function 


Sy(f)  «-  Ry(t1,t2)  (60) 

when 

Ry(t1,t2)  = R (t)  (61) 

Hard-Limiter . The  autocorrelation  function  of  the 
limiter  output  is  defined  as  the  statistical  average 

Ry(t1,t2)  = E(g[x(t1)].g[x(t2)]}  (62) 

where  t^  and  t2  are  two  arbitrary  times.  Since  the  ideal 
limiter  characteristic  as  described  by  Eq  (25)  can  be  repre- 
sented by  a sum  of  contour  integrals  in  the  complex  plane 
(w  = u + jv) , the  representation  is  given  by  (Ref  2: 277-282) 

g(x)  = yj7r{/c,  exp  (wx)  +/c_ex p (wx)  (63) 

The  autocorrelation  function  of  the  output  y of  a non-linear 
device  when  the  inputs  consist  of  sinusoidal  signals  plus 
independent  Gaussian  noise  is  given  by  (Ref  13:511): 


28 


r 


Ry(tl’t2>  ■ kf0Rn(tl-t2,^]y7  'c£<V  "2  ' 


exp  (o^  w^  / 2)  dw1/cf(w2)  w^  • 


.2  2 


exp  (a 2 w2  /2  dw2Mx(w1,w2) 


(64) 


where 

R 


n 


= autocorrelation  function  of  the  noise 


Mx (W1 * w2 ) = characteristic  function  of  the  input  signal 

f (w)  = bilateral  Laplace  transform  of  the  character- 

istic function  of  the  non-linear  device 


/ 


= integration  over  the  line  w = u - j<=°  to 
w = u + j°° 


The  normalized  transfer  function  for  a hard-limiter  output 
is  given  by  Eq  (25) . The  above  results  have  been  stated 
without  showing  the  actual  derivation.  A complete  analysis  of 
the  limiter  output  which  culminates  with  the  derivation  of  the 
autocorrelation  function  is  documented  in  (Ref  1) . 

After  considerable  manipulation  of  mathematical  equations 
and  analysis  of  those  equations,  Shaft  (Ref  5)  was  able  to 
derive  the  power  spectrum  by  transforming  the  autocorrelation 
function.  The  resulting  relationship  for  the  magnitude  of 
the  power  spectrum  at  the  various  frequencies  is  given  by 

Sml  ...mn=  {-°°  ^ Jml  (Alx)  ,,,Jmn(Anx) 


2 2 

exp  (■'■---x— ) dx) 


(65) 


r 


1 


largest  magnitude  are  those  which  satisfy  lm^lmax  2. 

The  subscripts  m.  . . .m  denote  the  harmonic  content  of  the 
in 

term  under  consideration.  For  example,  100... 0 indicates 
that  the  term  consists  of  the  fundamental  of  signal  (m  = 1) 
and  no  other  term  contributes  (m^  = 0,  i = 2...n).  Similarly! 
210... 0 is  the  intermodulation  term  consisting  of  the  second 
harmonic  of  signal  one  minus  the  fundamental  of  signal  two. 

The  condition 

n 

I m.  = 1 (66) 

i=l  1 


must  be  satisfied  for  a term  to  fall  in  the  first  zone. 

Soft-Limiter.  It  was  noted  by  Jones  (Ref  8),  that  if  a 
— 

soft-limiter  is  used,  the  exponent  in  Eq  (65)  must  be  replaced 
2 2 

by  [(a  + a ) x/2]  . The  variable  "a"  describes  the  slope  of 

the  error  function,  and  how  fast  saturation  is  reached.  The 
effect  of  adding  a soft-limiter  on  the  average  output  is, 
therefore,  mathematically  identical  to  adding  noise  to  a 
hard-limiter,  provided  that  the  error  function  representation 
can  be  used  for  the  limiter  characteristic. 

Approximation.  Since  Eq  (65)  is  very  difficult  to  evalu- 
ate numerically  for  large  n and  because  there  are  no  close 
form  solutions  for  n>2,  a useful  series  of  approximations  can 
be  derived  by  using  the  expression 


which  was  introduced  by  Gyi  and  used  by  Shaft  (Ref  14)  to 
produce  the  magnitude  of  the  desired  signals  and  unwanted 
signals . 

Two  Signals.  The  two  signals  under  consideration  are 
assumed  to  be  of  different  strength.  The  hard-limiter  output 
for  the  strong  term  is  expressed  as 

si o - (/:  va2x)  exp  dx}  (69) 

The  noise  in  this  case  was  assumed  to  be  negligible.  Hence, 
the  exponential  term  was  considered  to  be  unity.  After  several 
intermediate  steps  and  a change  of  variables,  y = A-^x  , and 
dy  = A-^dx  , the  output  is 

S10  = iji(y)  J0f^7)  dy}2  (69a) 

After  Eq  (14.4.33)  from  (Ref  1)  was  used  to  evaluate  Eq  (69a), 
the  results  obtained  were  manipulated  considerably  and  give 
the  approximate  term  for  the  strong  signal  output  as 

S10  as  1 (69b) 

The  strong  term  receives  most  of  the  power  from  the  non- 
linearity. An  approximate  expression  for  the  power  share  of 
the  weak  term  was  derived  using  the  same  technique  as  before. 
The  weak  term  output  is  approximately 

S01  “ T<jq>2  S10  <70> 

Its  power  level  is  6 dB  below  that  of  the  strong  term.  The 
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power  output  of  a hard-limiter  is  divided  unequally  among 
the  inputs,  with  the  stronger  input  being  allocated  a larger 
share  of  the  power  at  the  expense  of  the  weaker  signal. 

The  general  equation  for  the  output  from  a hard-limiter 
for  two  inputs  is  given  by  (Ref  14:658) 

/ (A2/A1)q  [ (r  + q)/2]  ^ 2 
Srq  = l 2r(q  + l)'r[  (r-q+Z)/2J  » (71) 

From  Eq  (71),  the  intermodulation  distortion  for  r = 2 and 
q = 1 is  found  to  be 


S 


21 


(72) 


Therefore,  intermodulation  distortion  can  be  obtained  and 
controlled  by  the  proper  selection  of  carrier  frequencies. 

The  center  frequencies  of  the  carriers  are  selected  in  a 
manner  that  insures  that  their  intermodulation  products  fall 
at  frequencies  out  of  the  desired  zone  or  where  they  minimally 
distort  the  desired  output. 

Several  Signals.  When  there  are  several  signals  present 
at  the  input  and  no  one  dominates,  the  noise  term  which  was 
previously  considered  negligible  and  ignored  must  now  be  in- 
cluded. The  output  term  for  one  of  the  carriers  is  given  by 

S100...0  - «'o  jFVV>  <--4  72^)2  C73J 


where 


= o2  ♦ 


? 1 A2 

2 7 Ai 

i = 2 1 1 


2 

°2 
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(74) 


2 

The  term  a ^ represents  noise  and  those  signals  whose  terms 
are  not  of  interest  are  also  considered  noise  and  are  summed 
into  the  noise  term  as  given  in  Eq  (74).  After  applying 
Eq  (11.4.28)  from  (Ref  1)  the  calculated  magnitude  for  the 
fundamental  frequency  of  the  first  signal  is  given  as 


0.7854 


(75) 


which  yields  a suppression  factor  of  -1.05  dB . This  is  a 
well  known  and  well  documented  result.  The  magnitude  for  the 
other  terms  at  their  fundamental  frequencies  is  also  given  by 
Eq  (73),  since  each  term  of  the  n terms  was  considered  to  be 
equal  in  magnitude  to  each  individual  term  in  the  sum. 

Using  the  approximation  introduced  and  defined  earlier  in 
Eq  (67),  the  desired  and  distortion  magnitude  terms  can  be 
obtained.  Those  calculated  results  allow  the  system  designer 
the  choice  to  select  appropriate  frequencies  to  eliminate 
distortion  terms  from  the  fundamental  zone  or  minimize  their 
effect.  The  total  intermodulation  power  cannot  be  altered, 
but  it  is  possible  to  redistribute  it  to  frequencies  where  it 
does  no  greater  harm. 

Intermodulation  Levels.  The  crossproducts  that  are  pro- 
duced when  several  signals  are  passed  through  a hard-limiter 
are  of  interest  since  they  represent  interfering  terms  or 
noise.  Values  for  the  magnitudes  have  been  obtained  from 
Eq  (65)  by  Shaft  (Ref  13).  Shaft  presented  his  numerical 
results  in  graphical  form.  His  results  were  extracted  and  are 
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Fig  9.  Magnitude  of  Individual  Terms 
Versus  the  Number  of  Signals  -n 


included  here  in  Fig  9.  The  magnitudes  shown  are  a function 
of  n signals,  all  of  which  have  equal  amplitude.  Several 
important  observations  can  be  noted  by  careful  inspection  of 
the  curves  in  the  previously  mentioned  figure.  These  obser- 
vations are  as  follow:  (1)  power  in  an  individual  term 
decreases  linearly  with  n;  (2)  total  signal  power  is  approxi- 
mately one  dB  less  than  the  total  power  and  remains  relatively 
constant  with  n;  (3)  power  in  an  individual  term  of  third  and 

fifth  order  decrease  as  n and  n , respectively;  and  (4)  band- 

% 

width  occupancy  for  the  signals  increases  linearly  with  n. 

Third  order  terms  are  defined  as 

N 


where  the  nus  are  integers  which  take  positive  and  negative 
values.  Antipodal  biphase-modulated  signals  (BPSK)  generate 
third  order  (2,1)  intermodulation  (IM)  products  which  are  also 
BPSK,  but  whose  frequency  is  offset.  Only  the  dominant  inter- 
modulation products  are  of  interest  because  they  occur  most 
frequently  and  have  four  times  the  power  of  the  210... 0 type 
IM  product. 

Number  of  Intermodulation  Products.  The  number  of  inter- 
modulation products  of  various  orders  are  tabulated  in  Table  I 
(Ref  15:243).  An  approximate  analysis  was  accomplished  on  a 
large  number  of  sinusoids  to  derive  the  total  number  of 
intermodulation  products  for  specific  orders.  The  notation 
A,  B,  C,...  represents  any  of  the  input  frequencies.  The 
asterisks  indicate  the  dominant  form  of  the  third  and  fifth 
order  cross-products. 


TABLE  I 

Number  of  IM  Products  of  Various  Orders 


IM 

Form 

Order 

m 

Number  of 
Frequencies 

Total  Number 
of  IM  Products 

2A 

- B 

3 

2 

N(N-l) 

*A 

+ B 

- C 

3 

3 

N(N-l) (N- 2) / 2 

3A 

- B 

- C 

5 

3 

N(N-l) (N - 2 ) / 2 

*A 

+ B 

+ C - D - E 

5 

5 

N(N-l) (N-2)  * 

(N-3) (N-4) / 1 2 

IV.  Detection  of  PSK  Signals 


Signals  passing  through  a linear  communications  channel 
are  corrupted  by  noise,  but  their  detection  is  much  easier  as 
compared  to  signals  which  are  passed  through  a non-linear 
system  that  introduces  still  more  noise.  The  equations  that 
relate  the  probability  of  error  for  making  a correct  decision 
for  both  a linear  and  a non-linear  device  are  provided  in 
this  section.  In  addition,  the  probability  of  error  for  a 
hard-limiter  and  a soft-limiter  are  found  to  be  expressible 
by  the  same  equation.  However,  the  softness  parameter  of  the 
soft-limiter  must  be  taken  into  account. 

Linear  PSK  Signal 

The  problem  of  communicating  one  of  a set  of  M specified 
signals  (s^(t)}  over  a channel  which  is  disturbed  by  additive 
white  Gaussian  noise  always  reduces  to  a corresponding  vector 
communication  problem.  The  transmitter  signals  are  represented 
by  M points  {s^}  in  an  N-dimensional  space.  The  optimum 
receiver  divides  the  signal  space  into  a set  of  M disjoint 
decision  regions  H^.  At  any  point,  q is  assigned  to  if, 
and  only  if,  (Ref  18:245) 

I q - s k | 2 - Nq  lnP[mk]  < | q-s± | 2 - NQ  lnP[m.] 

for  all  i t k (77) 

The  signals,  as  represented  in  Fig  10,  are  a binary 
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H 

l 

f-ig  10.  Decision  Space  for  Antipodal  Detection 

antipodal  set.  To  make  a decision,  the  decision  space  is 
divided  into  two  parts  by  a line  perpendicular  to  s^  - Sq. 

A complete  derivation  for  binary  and  M-ary  PSK  detection  is 
found  in  Van  Tress  (Ref  17:254  262),  and  Wozencraft  and 
Jacobs  (Ref  18:245-264).  The  error  probability  for  antipodal 
binary  signals  as  given  by  the  above  references  is 

Pe  = erfc[(Eb/No)1/2]  (78) 

Eq  (78)  was  numerically  evaluated.  Figure  11  shows  those 
numerical  results  in  graphical  form. 

Since  s^(t)  and  s2(t)  have  finite  energy  in  time  interval 
-T,  their  respective  energy  can  be  obtained  from 

i = 1,  2 (79) 


Fig  11.  Probability  of  Error  Versus 
Bit-Energy-to-Noise  Density 

2 

for  antipodal  signals,  E^  = E^  = A T/2.  The  energy  contained 
in  each  signal  is  identically  to  the  other.  For  a specific 
probability  of  error,  the  corresponding  signal-to-noise 
density  ratio  can  be  obtained  from  Fig  11  or  from  Eq  (78) . 

Limited  PSK  Signal 

A binary  phased-shif t-keyed  (BPSK)  signal  is  assumed  to 
get  from  the  transmitter  to  the  receiver  via  a limiting 
repeater  in  the  transmission  channel.  The  model  for  the 
transmission  channel  is  shown  in  Fig  12.  Additive  white 
Gaussian  noise  is  introduced  on  both  the  uplink  (repeater 
noise)  and  the  downlink  (receiver  noise). 

Hard-Limited.  The  input  of  the  limiter  in  Fig  12  consists 
of  a binary  phase-modulated  signal  of  frequency  wc  and  amplitude 
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Block  Diagram  of  PSK  System 


"A"  plus  zero-mean  stationary  Gaussian  noise  of  rms  value  a. 
Two  assumptions  are  made  at  this  point.  They  are  as  follows: 
the  noise  bandwidth  is  much  smaller  than  the  center  frequency 
of  the  filter,  and  the  bandpass  filter  preceding  the  limiter 
is  wide  enough  to  pass  the  signal  with  negligible  distortion. 
The  limiter  input  is  as  expressed  in  Eq  (27)  . During  any  bit 
interval,  the  information  bearing  parameter  takes  values  of 
either  0 or  tt , depending  on  whether  a zero  or  a one  is  being 
sent.  The  bandpass  limiter  is  assumed  to  be  ideal  in  the 
sense  that  its  fundamental-zone  output  given  by  Eq  (31). 

The  envelope  variations  are  removed  by  the  limiter,  but 
the  phase  (which  includes  the  information  signal)  remains 
undistorted.  The  signal  is  then  transmitted  to  the  receiver 
where  noise  is  added  to  it  on  the  downlink  at  the  receiver 
front  end.  Again,  the  downlink  noise  is  assumed  to  be  white 
Gaussian  noise  and  independent  of  the  uplink  noise.  Hence, 
the  receiver  input  may  be  expressed  as 

y^t)  = a cos(w(.t+e1)  + u(t)  cos(w(;t  + 01) 

- v (t)  sin(w(.t  + e1)  (80) 

where  the  noise  is  now  referenced  to  the  phase  of  the  informa- 
tion being  transmitted.  After  the  signal  is  multiplied  by 
2 cos(wct)  and  low-pass  filtered,  the  baseband  output  at  the 
sampling  instant  is  given  by 

z(t)  = k(t)  cos  (9  2 + 9 2^  (81) 
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where 


e2  " tan  -1t  aVitu(t)]  <«2> 

and 

k(t)  = {[a  + u(t)]2  + v2(t)}1/2  (83) 

The  amplitude  of  k(t)  of  the  signal  is  determined  by  the  ampli- 
fication at  the  repeater  and  by  the  downlink  losses.  The 
probability  density  function  of  the  sum  of  two  statistically 
independent,  identically  distributed  angles,  in  this  case 
0^  and  0 2»  each  representing  the  phase  of  a sinusoidal  signal 
plus  white  Gaussian  noise,  is  expressible  as  (Ref  6:37) 

m 

f(e  ) - * 1 x re?  ♦ i)  • 

m = l 

^F^(^-,  m+1,  -p2)  cos  m0^  for  i=l,2  (84) 

where  is  the  signal -to-noise  ratio  and  is  the  con- 

fluent hypergeometric  function. 

Error  Probability  Evaluation.  Since  the  receiver  output 
can  be  expressed  at  the  sampling  instant  as  the  sum  of  two 
statistically  independent,  identically  distributed  angles  as 
in  Eq  (81),  the  probability  of  error  can  be  determined  from 
the  characteristic  function  and  its  Fourier  transform.  The 
probability  density  function  of  the  two  independent  angles  is 
assumed  to  be  of  the  form  f ( 0 ) = f(-0)  . Also,  the  probability 
of  a one  or  zero,  which  can  be  represented  by  different 
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voltage  levels,  is  assumed  equally  probable.  Each  transmitted 
symbol  is  statistically  independent  of  its  predecessors.  Thus, 
it  will  be  sufficient  to  determine  the  error  probability  of 
a one  or  zero  since  they  are  equally  probable  and  even  sym- 
metric with  respect  to  0.  The  probability  of  error  in 
detecting  a one  or  zero  representation  is  equal  to  the  prob- 
ability that  the  receiver  output  is;  let  us  arbitrarily 
select  negative  at  the  sampling  instant.  Hence, 

Pe  = l+JzM  d*  (85) 

where  fz(z)  is  the  probability  density  function  of  the 
receiver  output  z. 

To  calculate  P , f (z)  must  be  determined  by  first  cal- 

C Z. 

culating  the  characteristic  function  of  z and  then  taking 
its  Fourier  transform.  Also,  the  Pg  can  be  obtained  directly 
from  the  characteristic  function.  The  characteristic  function 
is  defined  as 

Cz(v)  A E[e'ivz]  = fz(z)  ejvz  dz  (86) 

from  which  the  density  function  can  be  obtained  as 

fz(z)  = 7*  C Cz(v)  e’jVZ  dz  (87) 

If  z is  any  random  variable  of  probability  density  function 
f (z)  and  characteristic  function  C (v)  , its  cumulative  distri- 
bution function  P(x)  is  (Ref  6:37) 
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P(x)  = ll  f(z)  dz  = \ - I /*  Im[C2(v)  e'ivx]dv/v  (88) 

The  probability  error  was  found  to  be  independent  of  k(t) 
because  the  decision  is  only  made  on  the  sign  of  z.  After 
considerable  manipulations,  the  probability  of  error  as  a 
function  of  signal-to-noise  ratios  at  the  input  to  the  limiter 
and  at  the  input  to  the  receiver  was  derived  by  Jain  (Ref  7: 
627-628)  as 


where 


and 


Pe  = £[1  - (1-k2) 1/2  • Ie(k,x)] 


(89) 


(90) 


(91) 


The  Rice  function  written  as  an  integral  is 

Ie(k,x)  = fo  e_t  Vkt)  dt 


Alternately,  the  above  equation  can  be  represented  as  a 
series  as 

I (k,x)  = Z (k/ 2) 2n  (2n)l  An/n!  n! 
e n=0  n 


(92) 


where 


A « 1 - [1  + x + Z x2n/2n!]e*x 

n=0 


(93) 
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Figure  13  was  plotted  from  the  numerical  results  which 
were  obtained  from  evaluating  Eq  (92) . The  probability  of 
error  is  plotted  versus  the  signal- to-noise  ratios  at  the  limi- 
ter and  receiver  input.  The  numerical  results  obtained  from 
Eq  (92)  do  not  agree  with  the  results  obtained  by  Jain  (Ref  7: 
627).  The  two  graphs  do  resemble  each  other  to  a large  extent. 
Their  disparity  is  most  obvious  at  the  lower  values  of  error 
probability.  It  appears  that  the  graph  of  Fig  13  gives  values 
of  one  error  for  every  two  predicted  by  Jain's  graphical 
results.  Further  study  of  this  problem  was  not  carried  out 
because  it  was  determined  that  the  discovery  as  to  the  prob- 
able cause  of  the  disparity  would  not  provide  significant 
insight  into  the  actual  problem  being  investigated  in  this 
paper . 

Soft-Limited.  Figure  12  will  once  again  be  used  as  the 
system  model,  but  now  the  limiter  is  soft.  A soft-limiter  and 
its  characteristics  have  been  previously  defined  by  Fig  6. 

The  analysis  of  a soft-limiter  parallels  that  of  a hard-limiter. 
The  results  obtained  are  similar  and  are  summarized  in  this 
section. 

After  the  signal  is  multiplied  by  the  correctly  phased 
reference  signal  and  low-pass  filtered  to  remove  the  double 
frequency  components,  the  filter  output  is 

z(t)  = d cos^  + 02)  (94) 

where 

d = {[vx(a)  + u(t)]2  + v2(t)}1/2  (95) 
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The  amplitude  of  "d"  of  the  received  signal  is  additionally 
altered  by  the  amplification  in  the  repeater  and  the  downlink 
losses  after  the  repeater. 

Error  Probability  Evaluation.  Since  the  receiver  output 
can  be  expressed  at  the  sampling  instant  as  the  sum  of  two 
statistically  independent,  identically  distributed  angles  as 
in  Eq  (81),  the  probability  of  error  can  be  determined  from 
Eq  (89).  It  is,  therefore,  obvious  that  the  hard-limited  and 
soft-limited  channel  are  both  characterized  by  the  same 
expression. 

Once  the  signal - to-noise  ratio  has  been  determined,  the 
curves  which  were  calculated  for  the  hard-limiter  for  different 
probabilities  of  error  are  also  used  for  the  soft-limiter.  At 
small  SNRs,  the  curves  shown  in  Fig  2 of  (Ref  7:627)  indicate 
that  the  hard-limited  channel  yields  error  probabilities  which 
exceed  those  of  a linear  or  soft-limited  channel.  However,  a 
very  small  variation  in  the  signal-to-noise  ratio  in  that 
region  results  in  a drastic  change  in  the  P . For  large  SNRs, 
the  advantage  of  the  hard-limiter  no  longer  exists. 


V.  Link  Calculations 


A communications  link,  which  employs  dedicated  satellite 
channels  of  a multi-channel  transponder  for  data  transfer,  is 
limited  in  transmission  capability  by  one  of  several  constrain- 
ing factors,  one  of  which  usually  dominates.  These  factors  are 
as  follows:  (1)  earth  station  uplink  power;  (2)  satellite 
downlink  power;  (3)  satellite  and  earth  terminal  noise  levels; 
(4)  channel  bandwidth.  The  major  dominant  constraining  factors 
are  most  often  the  channel  bandwidth  or  downlink  SNR.  The 
equations  which  make  use  of  the  above  mentioned  parameters 
characterize  a communications  system. 


Communications  Equations 

The  design  of  any  telecommunications  system  makes  use  of 
the  communications  equations  which  are  stated  below  (Ref  10: 
8-15)  . 


ITT*  Wt  d)2  LoLpGr  rr 


S - (1  - L ) rec 


FT 

o 


V 1 r 


c _ r rec 

FT  Lm  1 T" 


where 

Lq  = Antenna  Pointing  Loss 

L = Polarization  Loss 
P 


(96) 

(97) 


(98) 
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Pt  = Total  Transmitted  Power 

L = Modulation  Loss 
m 

G^.  = Transmitter  Antenna  Gain 

G„  = Receiver  Antenna  Gain 
r 

\ = Wavelength  of  r-f  Carrier 

At  = Antenna  Capture  Area 

Most  often,  the  noise  Nq  is  referenced  to  the  input  stage  of 
the  receiver.  The  noise  as  referenced  to  the  input  at  the 
satellite  repeater  and  the  earth  terminal  receiver  are  quite 
different  values.  Usually  the  noise  at  the  earth  terminal 
receiver  is  far  less  than  at  the  satellite  receiver.  The 
parameters  K and  T are  defined  as 

K = 1.38  x 10"23  W/Hz  - °K  (99) 

known  as  Boltzmann’s  constant,  while  T is  given  in  degrees 
Kelvin.  The  remaining  parameters  used  in  the  communications 
equations  are  defined  as  follows: 

S/N  = Received  signal-to-noise  spectral  density  ratio 
in  the  data  channel 

Pc/Nq  = Carrier  signal-to-noise  spectral  density 
(X/4ird)2  = Space  loss  ; where  d « 3.563  x 10^  Km 

Performance  Parameters 

The  block  diagram  of  a PSK  system  which,  in  this  case, 
is  considered  to  be  communications  satellite  path,  is  shown  in 
Chapter  IV,  Fig  12.  The  transmission  begins  from  an  earth 
station,  is  radiated  through  space  to  the  satellite  repeater. 
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is  amplified  and  translated  in  frequency  at  the  repeater, 
and  is  finally  re-radiated  through  space  to  another  earth 
station.  Noise  corrupts  the  signal  on  the  uplink  and  downlink 
paths.  Not  only  does  the  transmission  suffer  corruption  from 
noise  sources  and  path  losses,  but  it  must  also  compete  with 
other  transmissions  for  a share  of  the  total  satellite  power. 
As  has  been  stated  before,  the  stronger  signals  receive  a 
larger  share  of  the  total  satellite  power  at  the  expense  of 
the  weaker  signals.  The  disparity  can  be  as  high  as  6 dB. 

The  critical  parameter  to  be  calculated  is  the  ratio  of 
received  carrier  power- to-noise  density  C/NQ.  The  required 
C/Nq  for  a given  bit  rate  R is  related  to  the  Eb/NQ  (energy- 
per-bit- to-noise  density)  by  the  expression  (Ref  15:176) 
r E. 

(jrO  = xr  + R + M dB-Hz  (100) 

o required  o 

where  R is  the  bit  rate  of  the  information  being  transferred, 
and  M is  almost  always  fixed  at  6 dB  as  a link  margin  for 
X-band  operation. 

Uplink.  The  signal  power  at  the  satellite  input  is  given 
by 

Ps  = ptGtGs(A/4TTd)2  (101) 

where 

Pt  = Earth  station  transmitter  power 
Gt  = Earth  station  Antenna  Gain 
Gg  = Satellite  Receiver  Antenna  Gain 
2 

The  factor  (A/4ird)  is  commonly  referred  to  as  the  free  space 
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loss.  The  antenna  gain  expression  is 


G = n4TT-^-  (102) 

n 

where  n denotes  an  arbitrary  antenna  and  An  an  arbitrary 
aperture.  The  efficiency  n is  normally  .55  for  a parabolic 
reflector  operating  in  the  X-band.  The  noise  level  at  the 
satellite  input  is  expressed  as 

Nq  = KT  (103) 

Thus,  the  carrier  to  noise  spectral  density  is  of  the  form 

ps  2 

FT  = PtGtGs(X/4lld)  /KT  (104) 

o 

A complete  uplink  relationship  is  stated  for  the  signal  to 
noise  spectral  density  by  Eq  (104). 

Downlink.  The  signals,  as  they  pass  through  the  non- 
linearity, are  altered  significantly  and  have  a different 
power-sharing  relationship  at  the  output.  Hence,  y is  intro- 
duced here  to  represent  the  suppression  or  enhancement  a 
particular  signal  experiences.  The  signal  power  at  the  earth 
station  receiver  input  is  obtained  from 

pr  = PtGtGs(X/4TTd)2  (105 

\ 

where  the  symbols  are  as  previously  defined  since  the  antennas 
used  to  receive  and  transmit  are,  in  most  instances,  one  and 
the  same  at  the  satellite  and  earth  station.  However,  Eq  (105) 
must  be  altered  to  include  a suppression  factor  and  some  amount 
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of  re-radiated  noise.  The  signal  and  noise  passing  through 
the  repeater  and  amplified  can  be  expressed  as 


Signal  Power  = yEIRP  (Ps/Ps+KT)  (106) 

Noise  Power  = y * EIRP  (KT/PR+KT)  (107) 

where  y is  the  suppression  factor  and  EIRP  is  the  power  of 
the  repeater  output  amplified  by  the  antenna  gain.  Hence, 
the  receive  power  is 

pR  - Y-GsPtsPt(X/4Trd)2.(Ps/Ps  + KT)  (108) 

where  P^  is  the  radiated  power  from  the  satellite.  The 
noise  level  at  the  earth  station  receiver  consists  of  re- 
radiated noise  and  downlink  noise  and  is  given  by 

Nr  = Y*GsPts(X/47rd)2(KT/Ps  + KT)  + KT  • (109) 

where  T'  is  the  noise  temperature  referenced  to  the  receiver 
input.  Normally,  the  first  term  of  NR  is  very  small.  In  the 
same  manner  that  the  desired  signal  suffers  a path  loss,  so 
does  the  unwanted  re-radiated  signal.  NR  is  the  sum  of  two 
independent  white  Gaussian  noise  sources.  The  satellite  link 
is  said  to  be  limited  whenever  the  noise  dominates.  The  more 
common  situation  is  the  downlink-limited  channel  where  the 
satellite  EIRP  is  the  limiting  factor.  The  carrier-to-noise 
density  for  that  specific  case  is  approximated  as 

p 

FT  * PtsYPtGs(X/4ird)2/KT'  (U°) 

R 
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r ^ 

The  parameter  y takes  into  account  such  factors  as  power 
sharing  and  suppression.  In  effect,  y represents  the  share 
of  the  satellite  power  a particular  signal  receives  in  the 
competitive  environment  of  several  signals. 


VI.  Configuration  Analysis 


To  understand  how  a system  functions  with  an  input  of 
several  signals  which  are  not  all  of  the  same  amplitude  and 
information  carrying  capacity,  it  was  first  necessary  to 
analyze  the  performance  of  the  system  with  only  one  access. 
Then,  the  system  was  analyzed  with  several  inputs,  all  of  which 
were  identical  in  bandwidth  requirement  and  amplitude.  Those 
results  have  been  previously  derived  and  summarized  in  the 
preceding  chapters.  In  this  chapter,  those  results  are  com- 
bined in  an  effort  to  understand  the  degraded  performance  of 
the  system  when  several  signals  of  different  amplitudes  and 
information  bearing  capacity  interact  as  they  pass  through 
the  non-linear  repeater.  The  FDMA/PSK  access,  when  <J>  takes 
values  of  0 and  it,  are  found  not  to  approach  a Gaussian  dis- 
tribution and  found,  therefore,  not  to  appear  as  Gaussian 
noise,  regardless  of  how  many  signals  are  taken  into  the  sum. 
That  analysis  was  accomplished  to  establish  a lower  bound  on 
the  suppression  factor.  However,  the  result  turned  out  to  be 
the  upper  bound.  The  least  suppression  a signal  experiences 
in  the  presence  of  several  FDMA/BPSK  signals  is  1.05  dB.  In 
reality,  the  suppression  factor  will  always  be  greater  than 
that.  FDMA/PSK  signals  with  <p  varying  uniformly  between  -n 
and  it  do  approach  a Gaussian  distribution.  The  FDMA/PSK 
signals  suppress  the  TDMA  signal  by  at  least  1.05  dB  and  much 
more.  In  turn,  the  TDMA  signal  and  n-1  FDMA/BPSK  signals 


53 


suppress  the  n-th  FDMA  signal  by  at  least  1.05  dB.  The  sup- 
pression that  each  of  the  signal  experiences  is  related  to 
their  individual  power.  Most  signals  are  suppressed  by  up 
to  6 dB. 


Time  Division  Multiple  Access 

In  TDMA,  the  signaling  waveforms  are  orthogonal  in 


time  because  each  terminal  has  exclusive  use  of  the  satellite 


during  periodically  recurring  time  intervals  or  slots.  The 
slots  and  slot  assignments  are  not  considered  or  discussed 
here  since  it  has  been  assumed  that  those  problems  will  be 
solved  in  the  design  of  the  TDMA  equipment  and  subsequent  to 
this  study. 

Because  each  earth  terminal  transmits  and  receives 


information  to  and  from  the  satellite  in  separate  non-overlapping 
time  slots,  the  generation  of  intermodulation  products  in  the 
non-linear  device  is  avoided.  However,  a hard-limiter  does 
introduce  some  signal  suppression.  The  greatest  suppression 
is  -1.05  dB  for  very  small  SNRs  and  0 dB  for  large  SNRs,  as 
has  been  previously  calculated  and  as  shown  in  Fig  5. 

Moreover,  the  TWT  is  operated  at  maximum  output  to  take  full 
advantage  of  the  limited  and  precious  satellite  power.  The  com- 
plete end-to-end  communications  channel  for  a hard-limited  or 
soft-limited  TDMA/PSK  communications  channel  can  be  accomplished 
with  the  aid  of  Figs  5,  8,  and  11,  or  the  equations  which  were 
used  to  generate  the  graphs  for  those  figures.  It  was  determined 
that  the  limiting  factor  for  a single  TDMA  access  for  a specific 
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probability  of  error  was  the  satellite  power  and  not  the 
channel  bandwidth.  A link  calculation  was  accomplished  to 
illustrate  this  point  and  is  provided  in  tabular  form  in 
Table  II.  Also,  the  most  information  that  can  be  effectively 
transmitted  through  a 60  MHz  channel  was  found  to  be  30  Mb/s 
for  BPSK  and  60  Mb/s  for  QPSK.  Figure  14  shows  three 
accesses  at  the  operating  frequency  and  bandwidth  of  the 
DSCS  III  transponder  II.  The  accesses,  as  they  appear  in  the 
frequency  domain,  are  indistinguishable  from  each  other  except 
for  the  bandwidth  requirements.  Figure  14a  shows  two  accesses 
which  could  be  considered  to  be  FDMA  and  one  access  which 
could  be  considered  TDMA.  It  should  now  be  obvious  that  FDMA 
and  TDMA  with  phase  shift  modulation  have  the  same  frequency 
spectrum.  However,  the  TDMA  access  will  almost  always  have  a 
wider  frequency  spectrum.  Again,  looking  at  the  figure 
closely  reveals  that  the  three  accesses  have  identical  power. 
That  is,  if  the  area  under  their  respective  curves  is  somehow 
calculated  by  some  integration  scheme,  they  will  all  have  the 
same  value.  The  spectrum  with  the  wider  bandwidth  will,  of 
course,  yield  a lower  energy-per-bit  term  since  more  informa- 
tion is  being  transferred  by  that  access.  The  figure  just 
below  the  previously  discussed  one  shows  two  interesting 
phenomena.  (1)  The  energy-per-bit  for  each  of  the  three  signals 
is  now  identical.  The  individual  power  in  each  of  the  r-f 
carriers  was  adjusted  to  insure  that  each  bit  received  the 
same  power  ratio  as  for  the  access  with  less  information 
bearing  capacity.  (2)  The  spectral  density  of  the  largest 
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TABLE  II 


Link  Calculation  for  a Single  TDMA  Access 


Uplink  Calculations 


EIRP  (dBW) 

Gs  (dB) 

System  Losses  (dB) 

Free  Space  Loss  (dB) 

Received  Signal  Level  (dBW) 

Noise  Density  (dBW/Hz) 

Carrier  Power  to  Noise  Density  (dB-Hz) 

Downlink  Calculations 
EIRP  (dBW) 

Gt  (dB) 

Free  Space  Loss  (dB) 

Received  Signal  Level  (dBW) 

Noise  Density  (dBW/Hz) 

Carrier  Power  to  Noise  Density  (dB  - Hz) 
Carrier  Power  to  Noise  Density  (Required] 


43.7  to 


-129.0  to 
66.6  to 


39,  43, 


73.7 

29.0 

-3.0 

201.7 

-99.0 

-195.6 

96.6 


45 

43.7 

202.3 

113.6 

202.9 

89.3 


For  an  error  probability  Pg  = 10~  with  Binary 
phase  shift  keyed  - E^/NQ  = 9.6  dB 
C/KT  = R + Eb/NQ  + M 

Information  Rate  and  Bandwidth  Requirements 
R = 89.3  - 9.6  - 6.0  = 73.7 
R = 23.44  Mb/s 
Bandwidth  = 46.88  MHz 


Available  Transponder  Channel  Bandwidth  = 60  MHz 


Fig  14.  Multiple  Access  Representation  in  the 
Frequency  Domain,  (a)  Unequal  Eb/N0  for  each 
the  Three  Accesses,  (b)  Equal  £^/N0for  each 
of  the  Three  Accesses 
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access  overlaps  into  the  frequency  space  of  the  other  two 
accesses.  That  distortion  is  crosstalk  into  the  other 
accesses  and  results  in  degraded  performance. 

Simultaneous  Multiple  Accesses 

In  the  following  section,  an  attempt  was  made  to  prove 
that  a satellite  system  which  includes  one  TDMA  access  and 
several  FDMA  accesses  can,  in  most  cases,  be  modeled  as  a 
hard-limiter.  It  is  known  that  TDMA  bursts  occur  intermit- 
tently. That  is,  each  terminal  is  allocated  a time  slot  for 
transmission  by  some  control  method,  but  that  does  not  guar- 
antee that  some  slots  will  not  remain  vacant.  With  knowledge 
of  such  occurrences,  it  was  possible  to  analyze  the  output 
spectrum  of  the  BPL  with  the  TDMA  access  active  and  inactive. 
When  the  TDMA  access  is  inactive,  the  FDMA  accesses  have 
complete  use  of  the  satellite  power.  It  should  be  remembered 
that  the  TDMA  access  vacancies  or  inactive  slots  occur  very 
infrequently,  but  they  do  occur.  How  the  output  power  of  a 
BPL  is  shared  by  n accesses  of  equal  magnitude  at  the  input 
to  the  repeater  has  been  determined  in  previous  chapters. 

The  main  problem  under  consideration  was  the  analysis  of 
simultaneous  non-homogeneous  multiple  accesses  sharing  the 
limited  satellite  power.  The  nature  of  TDMA  communications 
required  that  the  access  drive  TWT  into  saturation.  It  is 
necessary  to  operate  at  saturation  to  maintain  a uniform  out- 
put because  of  the  many  users  that  contribute  data  for  trans- 
mission through  the  satellite  during  each  second.  The  users, 
by  some  control  method  which  is  not  of  concern  here,  are 


58 


assigned  and  removed  from  the  system  on  a random  basis, 
based  upon  user  priority  and  availability  of  slots  for  data 
transmission. 

System  Model.  The  system  model  for  several  accesses, 
which  includes  a TDMA  access  and  several  FDMA  accesses,  is 
shown  in  Fig  15.  The  n individual  signals  are  competing  for 
a share  of  the  satellite  power.  Their  interaction  through 
a piecewise  limiter  was  analyzed.  The  limiter  representation 
is  shown  in  Fig  16.  Whenever  the  signal  and  noise  components 
exceed  e,  the  output  of  the  limiter  can  be  considered  identi- 
cal to  that  of  a hard-limiter  as  given  by  Eq  (25).  It  is, 
therefore,  apparent  that  whenever  the  magnitude  does  not 
exceed  e,  the  limiter  will  be  operating  in  the  soft  region 
and  can  be  considered  a linear  amplifier  of  gain  K. 

Normally  a TWT  operates  as  a soft-limiter  when  its 
operating  point  is  backed-off  B dB.  B is  usually  adjusted 
to  3 dB  or  more  for  FDMA  operation.  The  result  is  lower  inter- 
modulation products,  but  the  drive  level  of  the  n signals  must 
be  carefully  controlled.  For  large  fixed  ground  terminals, 
accurate  uplink  power  control  is  feasible.  For  mobile  tacti- 
cal terminals,  it  may  prove  exceedingly  difficult,  if  not 
impossible,  to  keep  the  aggregate  drive  level  to  the  limiter 
at  the  proper  point. 

s 

Piecewise  limiter.  The  clipping  effect  of  the  linear 
piecewise  limiter  was  determined  by  Spilker  (Ref  15:232). 

When  the  input  A does  not  exceed  e,  the  output  is  B.^  = A. 

For  A > e,  the  peaks  of  the  sinusoid  drive  the  limiter  into 


It 
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Fig  15.  System  Model  for  Several  Accesses 


Fig  16.  Soft-Limiter  (Represented  by 
Piecewise  Limiter) 


saturation  and  appear  clipped  off  at  the  output.  The  output 
is,  therefore,  the  sum  of  the  contributions  from  the  limiting 
and  from  the  linear  regions.  Thus,  the  output  is  given  as 


Bx  = ~ {[d-(|)2]1/2  + | sin-1 


dll) 


The  magnitude  of  the  output  resembles  that  of  the  error  function 
limiter  which  is  graphically  shown  in  Fig  7.  Whenever  the 

6 


r 


signal  power  is  much  greater  than  that  of  the  noise,  the 
limiters  mentioned  above  operate  very  much  like  hard- limiters . 

Signal  Suppression.  The  suppression  factor  experienced 
by  one  sinusoidal  signal,  which  can  be  the  sum  of  n weak 
signals,  in  the  presence  of  a strong  sinusoid  has  been  found 
to  be  6 dB.  Also,  the  suppression  factor  experienced  by  the 
weak  sinusoid  decreases  to  1.05  dB  when  the  strong  signal  has 
a Gaussian  probability  density.  These  two  cases  are  the  most 
and  least  a sinusoid  can  be  suppressed,  respectively.  These 
results  have  been  derived  in  the  hard-limiter  suppression 
section  and  are  in  complete  agreement  with  Spilker's  results 
(Ref  15:226) . 

For  two  sinusoidal  inputs  with  envelopes  A(t)  = A and 
B(t)  = B,  the  weak  signal  B is  suppressed  by  6 dB.  The  weak 
signal  B is  suppressed  by  a factor  of  2 in  amplitude  and  a 
factor  of  4 in  power.  This  suppression  factor  is  valid  even 
when  the  weak  signal  is  composed  of  the  sum  of  a large  number 
of  sinusoids. 

For  a strong  Gaussian  interference  and  a weak  sinusoid, 
the  suppression  factor  is  0.7854  (Ref  15:226-227).  The  result, 
as  stated,  applies  to  the  case  of  n FDMA  accesses.  Each  and 
every  individual  access  plays  the  role  of  the  weak  signal  and 
is  suppressed  with  respect  to  the  total  output  power  by  the 
suppression  factor  -R. 

A configuration,  which  consists  of  several  FDMA  accesses 
and  one  TDMA  access,  was  considered  to  approximate  a Rician 
interference  which  consists  of  a Gaussian  noise  source  and  a 
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sinusoid.  The  FDMA  accesses  were  taken  to  represent  the 
Gaussian  noise;  the  TDMA  access  was  taken  to  represent  the 
sinusoid.  The  Rician  probability  density  is  given  by 
(Ref  15:229) 

f(V  = An  e‘(An+2r)/2  IQ[(2r)1/2  AJ  (112) 

where  r = Pc/Pn  which  is  the  power  ratios  of  the  fixed 
envelope  interference  sinusoid  to  the  Gaussian  noise.  The 
Rician  density  is  graphed  in  Fig  17a  for  various  SNRs . The 
factor  that  a sinusoid  experiences  in  the  presence  of  the 
noise  through  a bandpass  hard- limiter  can  be  computed  from 

R = J (1+r) [e‘r/2  I0(£)]2  (113) 

which  was  also  numerically  evaluated.  The  results  obtained 
are  presented  in  graphical  form  in  Fig  17b.  As  expected,  the 
suppression  factor  varies  between  a minimum  value  of  1.05  dB 
and  a maximum  value  of  6 dB.  As  stated  earlier,  a weak 
signal  in  the  presence  of  a white  Gaussian  noise  source 
suffers  the  least  suppression  while  it  suffers  the  most  sup- 
pression in  the  presence  of  a strong  sinusoid. 

Gaussian  Approximation.  It  was  necessary  at  this  point 
to  prove  that  the  sum  of  several  FDMA  accesses  with  binary 

phase  modulated  information  impressed  on  them  do  constitute 

\ 

Gaussian  noise.  To  do  so,  it  was  first  assumed  that  x(t,<£)  is 
a stochastic  process  which  can  be  represented  as 

x(t,<|>)  = Ai  cos(wcit  + <J>i)  i = l,...n  (114) 
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Fig  17.  Rician  Interference,  (a)  Probability 
Density  Function  vs.  SNR,  (b) 
Suppression  Factor  vs.  SNR 
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where  A is  the  peak  value  of  the  ith  access  and  <f>^  is  the 
information  bearing  variable.  The  information  variable,  <f>^, 
was  assumed  to  take  the  values,  0 and  it,  with  equal  probabil- 
ity. Hence,  the  first  order  probability  density  of  the 
process  was  found  to  be  as  represented  graphically  in  Fig  18. 
Hence,  E[x(t)]  = 0,  while  the  variance  of  x(t)  was  found  to 

2 2 

be  = | A cos  w^t | . It  is  now  possible  to  represent  the 

sum  of  the  random  process  as 

N N 

z(t)  = E x-(t)  = A E cos  (w  - 1 + <f>,)  (115) 

i=l  1 i=l  C1  1 

where  the  A's  are  equal  in  magnitude,  and  <Jk  ' s are  identically 
independently  distributed  random  variables  with  pdf  as  pre- 
viously defined. 

For  a fixed  t,  the  first  order  pdf  of  the  new  process  can 
be  found  by  invoking  the  Central  Limit  Theorem.  Therefore, 
the  process  z for  a fixed  t has  mean  and  variance  defined  as 

z = *(t)l£ixed  t - N(0,  A1  Z cos'  wcit)  (116) 

The  variance  was  found  from  the  sum  of  the  individual  vari- 
ances of  x^'s.  Because  x^'s  are  orthogonal  with  zero  mean, 
then  they  are  also  uncorrelated. 

Faulty  Approximation.  The  analysis  as  accomplished  in 
the  preceding  section  to  approximate  Gaussian  noise  with  n 
FDMA  accesses  was  erroneously  based  on  the  Central  Limit 
Theorem.  A closer  look  reveals  that  in  lieu  of  the  Central 
Limit  Theorem,  the  Chernoff  bound  should  have  been  applied. 
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Fig  18.  Probability  Density  of  the  Stochastic 

Process  -x 

Using  the  Chernoff  bound,  it  was  proved  in  (Ref  18:108-109) 
that  the  sum  of  n random  processes  which  represent  BPSK 
modulated  r-f  carriers  do  not  approach  a Gaussian  distribution 
regardless  of  how  many  terms  are  considered  into  the  sum.  A 
specific  case  was  considered  to  prove  the  aforementioned 
point.  N binary  random  variables  were  considered  in  which 
each  assumes  the  value  of  0 or  1 with  equal  probability. 

The  difference  between  the  values  for  the  approximation  and 
the  exact  expressions  was  found  to  be  expressed  by  an  expo- 
nential function  of  N.  The  fractional  error  was  found  to 
become  larger  as  N became  larger.  Therefore,  the  sum  of  n 
BPSK  modulated  r-f  carriers  do  not  become  Gaussian  regardless 
how  many  terms  are  considered  into  the  sum. 

Without  doubt,  it  has  been  proven  that  the  sum  of 
several  FDMA  accesses  will  never  appear  as  Gaussian  noise. 
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Hence,  the  FDMA  signals  will  experience  a suppression  factor 
which  approaches  6 dB . The  suppression  factor  experienced  by 
the  TDMA  signal  in  the  presence  of  several  FDMA  signals  can 
be  found  by  computing  the  ratio  - Pc/Pn.  A specific  suppres- 
sion value  is  obtained  from  Fig  17b  for  a specific  signal- 
to-noise  ratio. 

Non-homogeneous  Multiple  Accesses 

A method  to  utilize  more  of  the  unused  transponder 
channel  bandwidth  is  outlined  below.  This  method  is  an 
approximation,  but  it  does  serve  its  purpose  well  when  applied 
to  the  environment  of  several  accesses,  all  of  which  are 
homogenous  with  the  exception  of  one.  Figure  19  shows  how 
more  use  can  be  made  of  a channel  when  the  TDMA  access  does 
not  make  complete  use  of  the  channel  and  when  the  satellite 
power  can  support  more  accesses.  The  frequency  of  the  FDMA 
carrier  closest  to  the  TDMA  access  is  found  from  the  following 
equation 


f2  = fx  - 3Af  (117) 

where 

B = Transponder  Channel  Bandwidth 

f^  = Carrier  Frequency  of  TDMA  access 

f_  = Carrier  Frequency  of  FDMA  access  closest  to  the 
TDMA  access 

Af  = Phase  Modulation  rate  of  the  TDMA  access 

Approximately  one-fifth  of  the  bandwidth  can  be  used  for 
the  FDMA  accesses.  The  channel  configuration  must  be  as 
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Fig  19.  Frequency  Separation  to  Minimise 
Intermodulations  on  the  Desired 
Frequency  Spectrum 

depicted  in  Fig  19.  The  dominant  intermodulation  products 
fall  outside  the  desired  frequency  spectrum.  Hence,  the 
noise  terms  of  the  intermodulation  products  are  placed  at 
frequencies  where  they  do  no  greater  harm. 
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VII.  Conclusion 


The  essential  background  information  and  facts  for  the 
analysis  of  the  problem  are  included  in  this  report  and 
summarized  in  this  chapter.  A recommended  solution  is  pro- 
vided, although  the  sponsoring  organization  is  provided  with 
sufficient  facts  to  evaluate  other  solutions.  First,  the 
proposed  use  of  the  system  is  restated.  Basically,  Hq  AFCS 
desires  to  implement,  if  it  is  deemed  feasible,  the  use  of 
simultaneous  non-homogeneous  accesses  through  a common 
transponder  channel. 


Summary 

Technology  has  progressed  to  the  point  where  it  is  now 
feasible  to  design  TDMA  equipment  that  fully  utilizes  the 
entire  power  and  bandwidth  of  a single  transponder  channel 
such  as  that  of  the  DSCS  III  communications  satellite  channel 
II.  A link  calculation  was  presented  in  Chapter  VI  which 
proved  that  the  SHF  Ground  Mobile  Forces  Terminals,  even 
when  equipped  with  eight-foot  parabolic  dishes  in  a tactical 
environment,  can  make  complete  use  of  the  single  transponder 
channel.  However,  the  system  used  was  one  with  a probability 
of  error  of  10  5 and  with  BPSK  modulation.  Moreover,  in  a 
tactical  environment,  it  is  extremely  difficult,  if  not 
impossible,  to  control  the  aggregate  drive  level  to  the  input 
of  the  TIVT.  Hence,  a single  input  which  drives  the  TWT  into 
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saturation  is  more  appropriate.  For  a specific  P , Fig  13 
in  Chapter  IV  can  be  used  to  find  the  required  SNR  at  the 
input  to  both  the  satellite  channel  and  the  ground  re- 
ceiver when  a single  TDMA/PSK  user  accesses  the  channel. 

The  deviation  from  a single  access  to  several  non- 
homogeneous  accesses  is  extremely  costly.  Beyond  doubt,  it 
was  proved  in  Chapter  III  that  these  accesses  will  experience 
both  suppression  and  intermodulation  distortion.  How  much 
suppression  and  intermodulation  distortion  is  experienced 
by  each  of  the  individual  signals  can  be  calculated  from 


Eqs  (65)  and  (113).  A method  to  derive  the  approximate  mag- 
nitude of  the  desired  and  undesired  terms  is  also  included. 

A specific  quantity  can  be  easily  calculated  if  a specific 
case  is  given. 

Several  individuals  who  have  published  articles  in 
journals,  for  example,  Shaft  (Ref  14),  have  shown  that  the 
strong  signal  (TDMA)  should  be  assigned  a higher  carrier 
frequency  than  that  of  the  FDMA  accesses.  In  addition,  Shaft 
recommends,  without  specifying,  that  the  individual  TDMA 
access  should  be  separated  by  a large  frequency  from  the  FDMA 
signals.  The  separation  should  be  large  enough  so  that  the 
sum  and  difference  of  the  carrier  frequencies  fall  outside 
the  desired  bandwidth  of  the  transponder  channel.  The  appli- 
cation, as  described  above,  is  shown  graphically  in  Fig  19. 
However,  even  after  the  concept  of  Fig  19  is  used,  only  a 
small  fraction  of  the  transponder  channel  can  be  used.  The 
distortion  effects  outweigh  any  benefits  derived  from  using 
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FDMA  signals  in  the  presence  of  a TDMA  signal. 

Finally,  it  was  shown  in  Chapter  VI,  Fig  17,  that  each 
of  the  several  non-homogeneous  signals  will  experience  at 
the  minimum  1.05  dB  of  suppression.  Because  the  n FDMA 
signals  and  the  TDMA  signal  cannot  be  approximated  by  a 
Gaussian  distribution,  the  suppression  factor  for  each  of  the 
n weak  signals  approaches  6 dB.  The  stronger  TDMA  signal 
will  experience  the  suppression  as  found  by  using  the  results 
available  in  Fig  17.  The  suppression  of  the  stronger  TDMA 
signal  depends  on  its  power  to  the  ratio  of  the  sum  power  of 
the  n FDMA  signals.  But  it  will  always  be  greater  than  1.05  dB. 
As  Figs  11  and  13  show,  a decrease  in  the  signal-to-noise 
ratio  results  in  an  increase  in  the  probability  of  error. 

Recommended  Solution 

Because  technology  has  progressed  to  the  point  where  it 
is  now  feasible  to  design  TDMA  equipment  that  fully  utilizes 
the  entire  bandwidth  and  power  of  a single  transponder  channel, 
complete  and  sole  use  of  the  channel  by  one  TDMA  access  is 
recommended.  Suppression  and  intermodulation  distortion 
effects  are  eliminated.  The  DSCS  III  transponder  was  designed 
around  this  concept--exclusive  use  of  a transponder  channel 
by  several  homogeneous  accesses  or,  more  appropriately,  a 
single  TDMA  access. 
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